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Abstract 

Paralleling our recent computationally-intensive (quasi-Monte Carlo) work for the case N = 4 
(quant-ph/0308037), we undertake the task for N = 6 of computing to high numerical accuracy, the 
formulas of Sommers and Zyczkowski (quant-ph/0304041 ) for the (N 2 — l)-dimensional volume and 
(N 2 — 2)-dimensional hyperarea of the (separable and nonseparable) N x N density matrices, based 
on the Bures (minimal monotone) metric — and also their analogous formulas (quant-ph/0302197) 
for the (non-monotone) flat Hilbert-Schmidt metric. With the same seven billion well-distributed 
("low-discrepancy") sample points, we estimate the unknown volumes and hyperareas based on five 
additional (monotone) metrics of interest, including the Kubo-Mori and Wigner-Yanase. Further, 
we estimate all of these seven volume and seven hyperarea (unknown) quantities when restricted 
to the separable density matrices. The ratios of separable volumes (hyperareas) to separable plus 
nonseparable volumes (hyperareas) yield estimates of the separability probabilities of generically 
rank-six (rank- five) density matrices. The (rank-six) separability probabilities obtained based on 
the 35-dimensional volumes appear to be — independently of the metric (each of the seven inducing 
Haar measure) employed — twice as large as those (rank-five ones) based on the 34-dimensional 
hyperareas. (An additional estimate — 33.9982 — of the ratio of the rank-6 Hilbert-Schmidt 
separability probability to the rank-^ one is quite clearly close to integral too.) The doubling 
relationship also appears to hold S(4.15for the N = 4 case for the Hilbert-Schmidt metric, but not 
the others. We fit simple exact formulas to our estimates of the Hilbert-Schmidt separable volumes 
and hyperareas in both the N = 4 and N = 6 cases. 



PACS numbers: Valid PACS 03.65.Ud,03.67.-a, 02.60.Jh, 02.40.Ky, 02.40.Re 



I. INTRODUCTION 

In Part I of the influential paper , "Volume of the set of mixed entangled states", 
Zyczkowski, Sanpera, Horodecki and Lewenstein considered the "question of how many 
entangled or, respectively, separable states are there in the set of all quantum states" . They 
cited philosophical, practical and physical reasons for doing so. They gave a qualitative 
argument Q, sec. III.B] - - contrary to their initial supposition — that the measure of 
separable states could not be strictly zero. There has since been considerable work 
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141 ]. using various forms of measures, to determine/estimate the 



"volume of separable states" , as well as the volume of separable and nonseparable states 
[15! , and hence probabilities of separability. One somewhat surprising development has 
)een the (principally numerical) indication — two independent estimates being 0.137884 
and 0.138119 (sec. IVI C 21 below) that the volume of separable states itself can take on 
a notably elegant form, in particular, — l)/3 ~ 0.138071, for the case of qubit-qubit 
pairs endowed with the statistical distinguishability metric (four times the Bures metric). 
(However, there seems to be a paucity of ideas on how to formally prove or disprove such a 
conjecture.) The research reported below was undertaken initially with the specific purpose 
of finding whether a putative comparably elegant formula for the volume of separable qubit- 
qutrit pairs might exist. We will report below (sec. IVI D|) the obtaining of certain possible 
formulas that fit our numerical results well, but none of such striking simplicity (nor none 
that extends it, in any natural apparent fashion). But we also obtain some new-type results 
of substantial independent interest. 



In a recent highly comprehensive analysis [l5( (cf. 17. |), Sommers and Zyczkowski ob- 
tained "a fairly general expression for the Bures volume of the submanifold of the states of 
rank N — n of the set of complex {(3 = 2) or real {(3 = 1) iV x N density matrices 

qi p,Bur es ) = dn ?r (d " +1)/2 r(j/3/2)r[l + (2n + j-l)/3/2] 

T(K + l)/2) * =1 T[(n + j)/3/2}T[l + (n + j-l)(3/2Y {) 

where d n = (N — n)[l + (N + n — l)/3/2] — 1 represents the dimensionality of the manifold 
. . . for n = the last factor simply equals unity and (|TJ gives the Bures volume of the entire 
space of density matrices, equal to that of a do- dimensional hyper-hemisphere with radius 
1/2. In the case n = 1 we obtain the volume of the surface of this set, while for n = N — 1 we 
get the volume of the set of pure states . . . which for /3 — 1(2) gives correctly the volume of 
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21j | . In a broader context, 



the real (complex) projective space of dimensions N — 1" 
spaces of density matrices (p) has been widely studied [1 
it serves as the minimal monotone metric |22| . 

In Part II of , Zyczkowski put forth a certain proposition. It was that "the link 
between the purity of the mixed states and the probability of entanglement is not sensitive 
to the measure [on the space of N x iV density matrices] used". His assertion was based 
on comparisons between a unitary product measure and an orthogonal product measure for 
the (qubit-qubit) case N = 4 ^, Fig. 2b]. The participation ratio — 1/Tr(p 2 ) — was used 
as the measure of purity. 

II. SEPARABILITY-PROBABILITY RATIOS 

In this study, we present (sec. I VI CJ) numerical evidence — limited largely to the specific 
(qubit-qutrit) case N = 6 — for a somewhat related proposition (which appears to be 



possibly topological in nature 



It is that a certain "ratio of ratios'' 
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is equal to 2, independently of the measure used — where the possible measures (including 
the just-discussed Bures) is comprised of volume elements (all incorporating Haar measure 
as a factor) of certain metrics defined on the N x N density matrices. Here by 

g(2, metric) 

Timetric __ NA (Q\ 

^sep+nonsep — (2,metric) ' ^ ' 

is indicated the ratio of the hyperarea of the (N 2 — 2)-dimensional boundary of the (iV 2 — 1)- 
dimensional convex set (CV) of N x N density matrices to the total volume of Cn- Further, 

j-, (2, metric) 

Tjmetric _ N,l / *\ 

n sep = „(2,metric) V / 

is the same type of hyperarea- volume ratio, but now restricted to the (classical/non- 
quantum) subset of Cat composed of the separable states (which we designate using 
E rather than S). A simple algebraic rearrangement of quotients then reveals that Q metric 
(J2J) is also interpretable as the ratio 

y-,[metric,rank— N] 

^metric _ N /r\ 

y-)[metric,rank—(N—l)] v ' 
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of the separability probability of the totality of (generically rank-iV) states in C 



N 



y-, (2, metric) 

j-)[metric,rank—N] N,0 /«\ 

N = q (2,metric) V ' 
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to the separability probability 

i(2, metric) 



p[metric,rariK—(l\ — 1)1 iv,i /-,\ 

AT = Q (2,metric) v'/ 

of the (generically rank-(iV — 1)) states that lie on the boundary of Cat. 



III. METRICS OF INTEREST 

Let us apply the Zyczkowski-Sommers Bures formula (0) to the two cases that will be of 
specific interest in this study, N = 6, n = 0, (3 = 2 and N = 6, n = 1, (3 = 2 — that is, 
the Bures 35-dimensional volume and 34-dimensional hyperarea of the complex 6x6 density 
matrices. (It would, of course, also be of interest to study the real case, (3 = 1, though we 
have not undertaken any work in that direction.) We then have that 

_18 

n{2,Bures) _ " ^ - 27075 10~ 1? (8) 

6,0 " 12221326970165372387328000 ~ K } 



and 

Sf[ Bures) = — « 2.04457 • 10" 15 . (9) 

6,1 138339065763438059520000 \j 

Here, we are able (somewhat paralleling our recent work for the qubit-qubit case N = 4 |3j , 
but in a rather more systematic manner ab initio than there), through advanced numerical 
(quasi- Monte Carlo/quasi-random) methods, to reproduce both of these values (jBJ), ©, to 
a considerable accuracy. At the same time, we compute numerical values — it would seem 
reasonable to presume, at least initially, with roughly the same level of accuracy — of these 
two quantities, but for the replacement of the Bures metric by five other monotone metrics 



of interest. These are the Kubo-Mori 



3, y, y, q> (arithmetic) average |3j , Wigner- 



Yanase QHHQ. Grosse-Krattenthaler- Slater (GKS) Q (or "quasi-Bures" Q) and 
(geometric) average monotone metrics — the two "averages" being formed from the minimal 
(Bures) and maximal (Yuen-Lax j^) monotone metrics, following the suggested procedure 
in [36J, eq. (20)]. No proven formulas, such as (JTJ), are presently available for these other 
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various quantities, although our research in [3] had suggested that the Kubo-Mori volume 
of the N x N density matrices is expressible as 

Q (2,KM) _ r)N(N-l)/2c{2,Bures) nf .N 

which for our case of N = 6 would give 

Sff M) = 32768^ ures) . (11) 



In light o 
monotone 



the considerable attention rece ntly devoted to the (Riemannian, but non- 



]) Hilbert-Schmidt metric [la ll2l , including the availability of exact volume 



and hypersurface formulas [l6[, we include it in supplementary analyses too. Further, we 
estimate for all these seven (six monotone and one non-monotone) metrics the (unknown) 
35-dimensional volumes and 34-dimensional hyperareas restricted to the separable 2x3 and 
3x2 systems. Then, we can, obviously, by taking ratios of separable quantities to their 
separable plus nonseparable counterparts obtain "probabilities of separability" — ja tpp_ic 
which was first investigated in pj, and studied further, using the Bures metric, in 

IV. TWO FORMS OF PARTIAL TRANSPOSITION 



We will employ the convenient Peres-Horodecki necessary and sufficient positive partial 
transposition criterion for separability [sol . 4fl| — asserting that a4x4or6x6 density matrix 
is separable if and only if all the eigenvalues of its partial transpose are nonnegative. (In the 
jit] case, it simply suffices to test the determinant of the partial transpose for 



4x4 [qubit-qu 

nonnegativity 4l|,|42(.) But in the 6x6 case, we have the qualitative difference that partial 
transposes can be determined in (at least) two inequivalent ways, either by transposing 
in place, in the natural blockwise manner, the nine 2x2 submatrices or the four 3x3 
sob— H e, (20)]. ( Obvi OUsly , such a oonu^eoe, arises in a bipartite Sy ste m 
onlv if the dimensions of the two parts are unequal.) We will throughout this study — as 
in y - - at the expense of added computation, analyze results using both forms of partial 
transpose. 

It is our anticipation — although yet without a formal demonstration — that in the limit 
of large sample size, the two sets of (separable volume and separable hyperarea) results of 
interest here should converge to true common values. Now, the author must admit that 
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he initially thought that it made no difference at all in which of the two ways the partial 
transpose was taken, that is, a 6 x 6 density matrix would either pass or fail both tests. Also, 
this seems to be a common attitude in the quantum information community (as judged by a 
number of personal reactions [cf. Q, fn. 2]]). Therefore, we present below a specific example 
of a 6 x 6 density matrix (pi) that remains a density matrix if its four 3x3 blocks are 
transposed, but not its nine 2x2 blocks, since the latter result has a negative eigenvalue 
(-0.00129836). 
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K. Zyczkowski has pointed out that the question of whether a given state p is entangled 
or not depends crucially upon the decomposition of the composite Hilbert space Ha ® H B 



(cf. 43j,|44|). For instance, for the simplest 2x2 case, the maximally entangled Bell state 



becomes "separable", he points out, if one considers entanglement with respect to another 
division of the space, e. g. A 1 = {$+, $-}, B' = So, it should not be surprising, 

at least in retrospect, that some states are separable with respect to one form of partial 
transposition, and not the other. In the course of examining this issue, we found that if 
one starts with an arbitrary 6x6 matrix (M), and alternates the two forms of partial 
transpostion on it, after twelve (=2x6) iterations of this process, one arrives back at 
the original 6x6 matrix. So, in group-theoretic terms, if we denote the three-by-three 



operation by a% and the two-by-two operation by a 2 , we have idempotency, a\ = a\ = I and 



la 2 a 3 



(a3d2) 6 = /• Further, one can go from the partial transpose 03 (M) to the partial 



transpose a 2 (M) via the matrix corresponding to the permutation {1, 4, 2, 5, 3, 6}. 
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Further, we constructed the related density matrix (p 2 ) 



1 - 

9 
















1 


n 
U 


1 


U 


n 
U 


23 




—i 
41 








1 

1 

5 





+ 

24 ~ 


I 

38 














1 

7 


1 

10 


l 

21 


—I 
13 





—i 

23 


1 1 

24 38 


10 ~ 21 


1 

6 










1 

41 





1 

13 







79 
630 


V 












/ 



Now, if p 2 is partially transposed using its nine 2x2 blocks, it gives the identical matrix as 
when pi is partially transposed using four 3x3 blocks. But the six eigenvalues of pi, that 
is, {0.322635,0.222222,0.1721,0.149677,0.119158,0.0142076} are not the same as the six 
eigenvalues of p 2 , that is, {0.300489,0.222222,0.204982,0.168304,0.0992763,0.00472644}. 
So, there can be no unitary transformation taking p\ to p 2 . (The possibility that p\ and 
P2 might have the same total measure(s) attached to them, can not formally be ruled out, 
however.) 



V. RESEARCH DESIGN 



Our main analysis will take the form of a quasi-Monte Carlo (Tezuka-Faure [45J, l^fj) 
numerical integration over the 35-dimensional hypercube ([0, l] 35 ) and a 34-dimensional 
subhypercube of it. In doing so, we implement a parameterization of the 6x6 density 
matrices in terms of thirty Euler angles (parameterizing 6x6 unitary matrices) and five 
hyperspherical angles (parameterizing the six eigenvalues — constrained to sum toi) HQ- 
We hold a single one of the five hyperspherical angles fixed in the 34-dimensional analysis, so 
that one of the six eigenvalues is always zero — and the density matrix is generically of rank 
five. The parameters are linearly transformed so that they each lie in the unit interval [0,1] 
and, thus, collectively in the unit hypercube. The computations consumed approximately 
five months using six PowerMacs in parallel, each generating a different segment of the 
Tezuka-Faure sequence. 
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A. Silver mean (\/2 — 1) conjectures for N = 4 



We have previously pursued a similar numerical analysis in investigating the separable 
and nonseparable volumes and hyperareas of the 4x4 density matrices Q|. Highly accurate 
results (as gauged in terms of known Bures quantities — based on two billion points 
of a Tezuka-Faure ("low discrepancy") sequence lying in the 15-dimensional hypercube - 
led us to advance several strikingly simple conjectures. For example, it was indicated that 
the Kubo-Mori volume of separable and nonseparable states was exactly 64 = 2 6 times the 
known Bures volume. (The exponent 6 is expressible — in terms of our general conjecture 
(jmj), relating the Bures and Kubo-Mori volumes — as N(N — l)/2, with N = 4.) Most 
prominently, though, it was conjectured that the statistical distinguishability (SD) volume 
of separable states is and lOu^g in terms of (four times) the Kubo-Mori metric. Here, 



V2 



0.414214 is the "silver mean" 



identically four times the Bures metric 



50 



51 



(cf. 53]). The SD metric is 



54|. (Consequently, the SD 15-dimensional volume 
of the 4x4 complex density matrices is 2 15 times that of the Bures volume — given by 
formula (0) for N = A,n = 0, /3 = 2 - - thus equaling the volume of a 15-dimensional 
hyper-hemisphere with radius 1, rather than ~ as in the Bures case itself jl^].) 

Unfortunately, there appears to be little in the way of indications in the literature, as 
to how one might formally prove or disprove these conjectures - - "brute force" symbolic 
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sec. 



integration seeming to be well beyond present technical/conceptual capabilities (cf. 
5]. (Certainly, Sommers and Zyczkowski 15] did not directly employ symbolic integration 
methodologies in deriving the Bures volume, hyperarea,...for iV-level [separable and non- 
separable] systems, but rather, principally, used concepts of random matrix theory.) One 
approach we have considered in this regard |7j is to parameterize the 15-dimensional convex 
set of bipartite qubit states in terms of the weights used in the expansion of the state in 



some basis of sixteen extreme separable 4 x 4 density matrices (cf. 55j). For a certain basis 
composed of SU (4) generators |£f 



n 



. l58| . the associated 15 x 15 Bures metric tensor 
is diagonal in form (having all entries equal) at the fully mixed state p, sec. II. F]. (Also, we 
have speculated that perhaps there is some way of "bypassing" the formidable computation 
of the Bures metric tensor, and yet being able to arrive at the required volume element.) 
Perhaps, though, at least in the Bures/ minimal monotone case, a proof might be based on 
the concept of "minimal volume" 



s 



B. Formulas for monotone metrics 



The monotone metrics (of which we study five, in addition to the Bures) can all be 
expressed in the general form 

g p {X',X) = h: a ^\{a\X\f3)\ 2 c monotone {\ a ,\p) (14) 



(cf. [18, 19]). Here X, X' lie in the tangent space of all Hermitian N x N density matrices 
p and | a), a = 1, 2 . . . are the eigenvectors of p with eigenvalues X a . Now, c monatone (X a , Xa) 
represents the specific Morozova-Chentsov function for the monotone metric in question [36]. 
This function takes the form for: (1) the Bures metric, 

2 

CBures(X a , A/3) = ■ — —J (15) 

*a + A/3 

(2) the Kubo-Mori metric (which, up up to a scale factor, is the unique monotone Riemannian 



metric with respect to which the exponential and mixture connections are dual 



log A a - log A/3 

CKM{X a ,X(3) = ; (16) 

A a — As 



(3) the (arithmetic) average metric (first discussed in 



, _ 4(A a + A/3) 

c arith (X a , A,) - A 2 +6AqA/? + A 2. (17) 



(4) the Wigner-Yanase metric (which corresponds to a space of constant curvature 29]); 

4 

cn/y(A Q , A/3) = 77=— — 7^=7! (18) 

(v x a + ^Xpy 

(5) the GKS/quasi-Bures metric (which yields the asymptotic redundancy for universal 
quantum data compression [33]); 

Cgks(K,X/3) = — : e; (19) 

A/3 

and (6) the (geometric) average metric (apparently previously unanalyzed), 

Cgeom (Xa, A/3) = -==. (20) 

V X a Xp 

(The results obtained below for the geometric average monotone metric seem, in retrospect, 
to be of little interest, other than indicating — that like the maximal monotone (Yuen-Lax) 

n 

metric itself |3J — volumes and hyperareas appear to be simply infinite in magnitude.) 
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FIG. 1: Ratios of the cumulative estimates of the 35-dimensional volume S, 



,(2,Bures) 
'6,0 



to its known 



value (jHj). For each addtional point shown — as in all the subsequent plots — ten million (10 7 ) 
values of the particular seven-billion-point Tezuka-Faure sequence have been generated. 

VI. ANALYSES 

A. Volumes and Hyperareas Based on Certain Monotone Metrics 

Using the first seven billion points of a Tezuka-Faure sequence, we obtained the results 
reported in Tables Ul lIXI and Figs. HllT^l We followed the Bures formulas in sees. III.C, 
III.D], substituting for (|T5j) the Morozova-Chentsov functions given above fT6|) - ([20|) to obtain 
the non-Bures counterparts. 

In Fig. Q we show the ratios of the cumulative estimates of the 35-dimensional volume 
S^Q Bures ' > to its known value (jHJ. Each successive point is based on ten million (10 7 ) more 
systematically sampled values in the 35-dimensional hypercube than the previous point in 
the computational sequence. In Fig. El we show the ratios of the cumulative estimates of the 
34-dimensional hyperarea sf^ ures ^ to its known value ©• Each successive point is based on 
ten million more sampled values in the 34-dimensional hypercube than the previous point 
in the computational sequence. The single Tezuka-Faure sequence we employ for all our 
purposes, however, is specifically designed as a 35-dimensional one — of which we take an 
essentially arbitrary 34-dimensional projection. This is arguably a sub o ptim al strategy for 
generating well-distributed points in the 34-dimensional hypercube (cf. [62, sec. 7]), but it 
is certainly highly computationally convenient for us (since we avoid having to generate a 
totally new 34-dimensional sequence — which would, we believe, increase our computation 
time roughly 50 percent), and seems to perform rather well. (In fact, as discussed below, 
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FIG. 2: Ratios of the cumulative estimates of the 34-dimensional hyperarea s& Bures ^ ^ Q known 
value © 




FIG. 3: A pair of cumulative estimates of the 35-dimensional Bures volume of separable qubit-qutrit 
states based on the two distinct forms of partial transposition 

the bias of our estimates seems to be — contrary to expectations — markedly less for the 
known [Bures and Hilbert-Schmidt] 34-dimensional hyperareas than for the 35-dimensional 
volumes.) 

We also present a joint plot (Fig. OJ) of the two sets of cumulative estimates of the Bures 
volume of separable qubit-qutrit states based on both forms of partial transposition. The 
estimates obtained using the four blocks of 3 x 3 submatrices, in general, dominate those 
using nine blocks of 2 x 2 submatrices. 

In TablelU we scale the estimates (which we denote using S) of the volumes and hyperareas 
by the known values © , (jUJ) of Sf^^^ and s& Bures *> ; w hile in Table |Il] we scale these 
estimates by the estimated values (7.22904- 10 -17 and 2.03991 ■ 10 -15 ) of these two quantities. 
(We use both approaches because we are uncertain as to which may be more revealing as to 
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TABLE I: Scaled estimates based on the Tezuka-Faure sequence of 7 billion points of the 35- 
dimensional volumes and 34-dimensional hyperareas of the 6x6 density matrices, using several 
monotone metrics. The scaling factors are the known values of the volume and hyperarea for the 
Bures metric, given by and more specifically for the cases: N = 6; n = 0, 1; and (3 = 2 by © 
and Q. 

, . Q(2,metric) / q{2, Bures) Q(2,metric) i q(2, Bures) 

metric b eo /<-'6,o ^6,1 /'-'6,i 

Bures 0.996899 0.999022 

KM 32419.4 45.4577 

arith 621.714 31.291 

WY 131.711 9.76835 

GKS 12.4001 3.55929 

geom 2.80011 • 10 44 1.44011 • 10 14 



TABLE II: Scaled estimates based on the Tezuka-Faure sequence of 7 billion points of the 35- 
dimensional volumes and 34-dimensional hyperareas of the 6x6 density matrices, using several 
monotone metrics. The scaling factors are the estimated values (S^Q Bures ^ = 7.2482 • 10~ 17 and 
g^^ wes ) _ 2.04257 • 10~ 15 ) of the volume and hyperarea for the Bures metric. 

, . Q(2,metric) , q(2, Bures) 0(^1 metric) 1 q(2, Bures) 

metric d 60 /^6fl b 61 /<-'6,i 

KM 32520.3 45.5022 

arith 623.648 31.3216 

WY 132.121 9.77791 

GKS 12.4387 3.56278 

geom 2.80882 • 10 44 1.44152 • 10 14 



possible exact formulas — an approach suggested by our work in [3j].) The results for the 
geometric average monotone metric in Table H] appear to be divergent. We might speculate 
that the middle four scaled hyperareas in the last column of Table |J correspond to the actual 
values 7 • 13/2 = 45.5, 2 2 • 5 2 /3 w 31.333, 3 • 13/4 = 9.75 and 7/2 = 3.5, and for the second 
column that we have 132 = 12-11 and 12, as actual values. 

In Tables IIHI and Tables IIV1 we report our estimates (scaled by the values obtained for 
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the Bures metric) of the volumes and hyperareas of the 6x6 separable complex density 
matrices. Let us note, however, that to compute the hyperarea of the complete boundary 
of the separable states, one must also include those 6x6 density matrices of full rank, 
the partial transposes of which have a zero eigenvalue, with all other eigenvalues being 
nonnegative |63J. (We do not compute this additional contribution here — as we undertook 
to do in our lower- dimensional analysis Q] — as it would slow quite considerably the overall 
process in which we are engaged, since high-degree polynomials would need to be solved at 
each iteration.) 

In (sj, we had been led to conjecture that that part of the 14-dimensional boundary of 
separable 4x4 density matrices consisting generically of Tank-four density matrices had SD 
hyperarea 55 g 9 A9 and that part composed of rank-three density matrices, 3 ^ 3 , for a total 14- 
dimensional boundary SD hyperarea of 98 ^ Ag . We, then, sought to apply the "Levy-Gromov 
isoperimetric inequality" 64| to the relation between the known and estimated SD volumes 
and hyperareas of the separable and separable plus nonseparable states J^, sec. VII. C]. 



Restricting ourselves now to considering only the separable density matrices, for Table ITTT1 
we computed the partial transposes of the 6x6 density matrices by transposing in place the 
four 3 x 3 sub matrices, while in Table ITVl we transposed in place the nine 2x2 submatrices. 

In Table |Vj we only require the density matrix in question to pass either of the two 
tests, while in Table IVI1 we require it to pass both tests for separability. (Of the 7 billion 
points of the Tezuka-Faure 35-dimensional sequence so far generated, approximately 2.91 
percent yielded density matrices passing the test for Table I, 2.84 percent for Table II, 4 
percent for Table III and 1.75 percent for Table IV. K. Zyczkowski commented that "it is not 
reasonable to ask about the probability that both partial transpositions are simultaneously 
positive, since one should not mix two different physical problems together".) 

In Table VII, we "pool" (average) the results for the separable volumes and hyperareas 
reported in Tables III and IV, based on the two distinct forms of partial transposition, to 
obtain possibly superior estimates of these quantities, which presumably are actually one 
and the same independent of the particular form of partial transposition. 

In Fig. H| we show the ratios of Sqq K to its conjectured value (fTTj) of 327685 , g 2 Q Bures ^. 
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TABLE III: Scaled estimates based on the Tezuka-Faure sequence of 7 billion points of the 35- 
dimensional volumes and 34-dimensional hyperareas of the separable 6x6 density matrices, using 
several monotone metrics. The scaling factors are the estimated values (Sg 2 Q Bnre ^ = 1.0739 • 10~ 19 
and f^'Bures) _ i_53g32 . i o -18 ) — the true values being unknown — of the separable volume and 
hyperarea for the Bures metric. To implement the Peres-Horodecki positive partial transposition 
criterion, we compute the partial transposes of the four 3x3 submatrices (blocks) of the density 



matrix. 


metric 


Bures-scaled separable volume 


Bures-scaled separable (rank-5) hyperarea 


KM 


8694.79 


9.43481 


arith 


220.75 


10.6415 


WY 


55.3839 


4.13924 


GKS 


7.97798 


2.28649 


geom 


3.33872 • 10 32 


3.61411 • 10 8 




FIG. 4: Ratios of the cumulative estimates of the 35-dimensional volume S^q K to its conjectured 
value of 327684 2 6 B " res) 

B. Volumes and Hyperareas based on the Hilbert-Schmidt Metric 

Along with the computations based on six distinct monotone metrics, reported above in 
sec. IVI Al we have at the same time carried out fully parallel analyses of the (Riemannian, 



but non-monotone) Hilbert-Schmidt metric 



371 ] . These have only been conducted after an 



earlier less-extensive form of this analysis [9 1 , reporting initial numerical estimates for the 
same six monotone metrics based on 600 million points of a Tezuka-Faure sequence, was 
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TABLE IV: Scaled estimates based on the Tezuka-Faure sequence of 7 billion points of the 35- 
dimensional volumes and 34-dimensional hyperareas of the separable 6x6 density matrices, using 
several monotone metrics. The scaling factors are the estimated values (Sg 2 Q BMres ' ) = 9.54508 • 1CT 20 
and yj( 2 ^ Bures ) _ 1.40208 • 10~ 18 ) — the true values being unknown — of the separable volume and 
hyperarea for the Bures metric. To implement the Peres-Horodecki positive partial transposition 
criterion, we compute the partial transposes of the nine 2x2 submatrices (blocks) of the density 
matrix. 



metric 



Bures-scaled separable volume 



Bures-scaled separable (rank-5) hyperarea 



KM 
arith 

WY 
GKS 
geom 



6465.86 
218.602 
55.5199 
7.92729 
5.4299 • 10 35 



9.0409 
10.3248 
4.05201 
2.26136 
4.35667 • 10 9 



TABLE V: Scaled estimates based on the Tezuka-Faure sequence of 7 billion points of the 35- 
dimensional volumes and 34-dimensional hyperareas of the separable 6x6 density matrices, using 
several monotone metrics. The scaling factors are the estimated values (1.99772- 10~ 19 and 2.90956- 
10~ 18 ) — the true values being unknown — for the Bures metric. A density matrix is included here 
if it passes either form of the positive partial transposition test. 



metric Bures-scaled separable volume Bures-scaled separable (rank-5) hyperarea 

KM 7735.7 9.30446 

arith 221.689 10.5467 

WY 55.8928 4.11453 

GKS 7.99075 2.28089 

geom 2.59619 • 10 35 2.29051 • 10 9 



posted. At that stage of our research, we had — with certainly some regrets — decided to 
fully redo the computations reported there. This was done to avoid a (somewhat inadvertent) 
programming limitation (which seemed of minor importance at the time) — a consequence 
essentially only of our, in time, having understood how to greatly speed up the computations 
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TABLE VI: Scaled estimates based on the Tezuka-Faure sequence of 7 billion points of the 35- 
dimensional volumes and 34-dimensional hyperareas of the separable 6x6 density matrices, using 
several monotone metrics. The scaling factors are the estimated values (3.06807- 10" 21 and 3.78991- 
10~ 32 ) — the true values being unknown — for the Bures metric. A density matrix is included 
here only if it passes both forms of the positive partial transposition test. 



metric 


Bures-scaled separable volume 


Bures-scaled separable 


(rank-5) hyperarea 


KM 
arith 

WY 
GKS 
geom 


1800.19 
92.7744 
26.4785 
5.56969 
3.96779 • 10 27 




3.99932 
5.3548 
2.55592 
1.77049 
1.09937 • 10 7 


TABLE VII: Scaled estimates obtained by pooling the results from Tables III and IV — based on the 
two forms of partial transposition — for the separable volumes and hyperareas. The Bures scaling 
factors (pooled volume and hyperarea) are Sg 2 Q Sures ^ = 1.0142- 10~ 19 and f^^Bures) _ 14707.PQ-I8 


metric 


Bures-scaled separable volume 


Bures-scaled separable 


(rank-5) hyperarea 


KM 
arith 

WY 
GKS 

geom 


7645.92 
219.739 
55.4479 
7.95413 
2.55692 • 10 35 




9.24704 
10.4905 
4.09766 
2.27537 
2.26584 • 10 9 



- of not being able to sample more than two billion Tezuka-Faure points. This fresh 
beginning (incorporating a much larger limitation, of which we here take advantage) allowed 
us then, as well, to additionally fully include the Hilbert-Schmidt metric. (It is somewhat 
unfortunate, however, at this point, that we had not conducted analyses based on the HS 
metric for the N = 4 qubit-qubit case, having restricted our earlier attention to monotone 
metrics only 3] [cf. sec. I VI C2\ .) 

Prior to Sommers and Zyczkowski reporting their exact formula (0) for the Bures volume 
of the submanifold of the states of rank N — n of the set of complex {(3 = 2) or real (/3 = 1) 
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N x N density matrices, they had obtained fully analogous formulas for the Hilbert-Schmidt 
metric, which for the specific volume (n = 0) case gives eq. (4.5)], 

= VNfryw-M * 11 ];-™ (2i) 



and hyperarea (n = 1) case Q, eq. (5.2)] gives, 

5 6>1 - viv i w r(iv)r(iV2 - 1) • 

For the (qubit-qutrit) case N = 6 under study in this paper, these give us, for the 
3 5- dimensional HS volume, 

1520749664069126407256340000000^ 
and for the 34-dimensional HS hyperarea 



15 

Sff S) = j= « 7.69334 10" 24 ( 23) 



15 

S ( ei HS) = = « 1.47483 10~ 21 . (24) 

8689998080395008041464800000 VE 

So, as above, using the Bures metric, we can further gauge the accuracy of the Tezuka-Faure 
numerical integration in terms of these known volumes and hyperareas. (This somewhat 
alleviates the shortcoming of the Tezuka-Faure procedure in not lending itself to statistical 
testing in any straightforward manner.) 

The estimated probability of separability is greater for the HS-metric than for any mono- 
tone one. (The minimal monotone or Bures metric appears to give the greatest probability 
in the nondenumerable class of monotone metrics. Also, the maximal monotone metric 
seems to give a zero probability 3].) Therefore, one might surmise that the much-discussed 
estimates of the sizes of the separable neighborhoods [n, 12, Q — which usually appear 
to be based on the HS or Frobenius metric — surrounding the fully mixed state are on the 
rather generous side (cf. |3]), relatively speaking. 

In Fig. Elwe show — paralleling Fig. [T] — the ratios of our cumulative estimates of S$ q 
to its known value (|23j). In Fig. El we show — paralleling Fig.|21 — the ratios of the cumulative 
estimates of S^{ HS ^ to its known value (|24p . 

A plot (Fig. Q) of the cumulative estimates of the Hilbert-Schmidt volume of separable 
qubit-qutrit states (for the two forms of partial transposition) is also presented. (The ratio of 
the two cumulative estimates at the final [seven billion] mark is 1.03236, while the comparable 
ratio is 1.12508 in the analogous Bures plot [FigJSJ.) 
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FIG. 5: Ratios of the cumulative estimates of the 35-dimensional Hilbert-Schmidt (Euclidean) 
volume Sj?Q HS ^ to its known value 
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FIG. 6: Ratios of the cumulative estimates of the 34-dimensional Hilbert-Schmidt hyperarea Sg { 
to its known value (I24[) 
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FIG. 7: Cumulative estimates of the 35-dimensional Hilbert-Schmidt volume of separable qubit- 
qutrit states (for the two possible forms of partial transposition) 
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In their two studies ^| , deriving exact formulas for the Bures and Hilbert-Schmidt 
volumes and hyperareas of the N x N density matrices, Sommers and Zyczkowski also 
explicitly derived expressions for the ratios of (iV 2 — 2)-dimensional hyperareas to (iV 2 — 1)- 
dimensional volumes. These were eq. (4.20)] 

s ( Nf ures) 2 r(iV 2 /2) 

lBures,N ~ ^Bures) ~ ^ T{N 2 /2 - 1/2) 1 ' 



and 



id eq. (6.5)] 

1hs,n = -gbj = v^(iV-l)(iV 2 - 1). (26) 

In the N = 6 Bures case, this ratio (equivalently what we have earlier denoted Rf^+nonsep) 
is 

S ?i BUTeS) 2 34 17179869184 

iBures 6 = ~TTB \" = -5 = ~ 28.1205, 27 

iBures,b {2 ,Bures) 3 2 - 5 ■ 11 • 19 • 23 • 29 • 31tt 194467185tt ' v ; 

and in the Hilbert-Schmidt case (equivalently R^ e p +nonsep ) , 

g(2,HS) 

insfi = -g^y = 35^ « 191.703. (28) 

^6,0 

(The Bures ratio grows proportionally with the dimensionality (D = N 2 — 1) of the N x N 
density matrices as D (for large N) jl^ . sec. IV. C], and as D 3 ^ 2 for the Hilbert-Schmidt 
ratio 16J, sec. VI].) Our sample estimates for these two quantities are 28.1804 and 192.468, 
respectively. In Table IV1111 we report these estimates, as well as the sample estimates for 
the other five metrics under study here. We also list the two known values and also give the 
corresponding ratios of hyperarea to volume for a 35-dimensional Euclidean ball having: (1) 
the same volume as for the metric in question; and (2) the same hyperarea. Only for the 



(flat) HS-metric are these last two ratios less than unity (cf. (l6j, sec. VI]). 

C. Separability-probability ratios 

1. The N = 6 qubit-qutrit case 

In Table HXl we list for the seven metrics the estimated ratios, which we denote i?™p iric , 
of the hyperarea (consisting of only the rank-five but not the rank-six 6x6 density matrices 
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TABLE VIII: Sample estimates of the ratio (Rf etric 



-,(2, metric) i q(2, metric) \ 



of the 34- 



L sep+nonsep 6,1 / 6,0 

dimensional hyperarea to the 35-dimensional volume for the seven metrics under study, and the 
corresponding ratios for a 35-dimensional Euclidean ball having: (1) the same volume as for the 
metric; and (2) the same hyperarea 

metric known ratio sample ratio (R^p%nonsep) isovolumetric ratio isohyperarea ratio 



Bures 
KM 
arith 
WY 
GKS 
geom 
HS 



28.1205 



191.703 



28.1804 
0.0394299 
1.41531 
2.08556 
8.07163 
1.44625 • 10~ 29 
192.468 



2.34553 
1245.79 
38.858 
27.5655 
7.61987 
2.45463 • 10 29 
0.543466 



2.40508 
1536.34 
43.2743 
30.39 
8.08886 
1.79638 • 10 30 
0.533806 



TABLE IX: Sample estimates of the ratio (i?™p*" c ) of the 34-dimensional hyperarea consisting 



only of rank-five 6x6 separable density matrices to the 35-dimensional separable volume for the 
seven metrics under study. In the last column there are given the ratios of ratios (O me *" c ) of the 
middle (third) column of Table IVIIII to these values 



metric 



nmetric 
^sep 



Q-metric — rpmetric J nmetric p 

J ' — ^sep+nonsepl ^sep r t 



metricfi] , j-^metricfi] 



Bures 
KM 
arith 
WY 
GKS 
geom 
HS 



14.501 
0.0175377 
0.692291 
1.07164 
4.14819 
1.28502 • 10~ 25 
94.9063 



1.94334 
2.24829 
2.04439 
1.94613 
1.94582 
0.000112547 
2.0279 



constituting the boundary of the separable density matrices j&l) to the volume of all the 
separable states themselves. We see that R S ^ Y ls quite close to 1. (The Wigner-Yanase 
metric is one of constant curvature js^.) In Fig.|H]we show the deviations of the cumulative 
estimates of R^ Y fr° m 1- 111 the last column of Table llXl there are given the ratios of ra- 
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FIG. 8: Deviations of the cumulative estimates of i?wy from 1. 

tios [} metrtc = R™J^+n 0nsep / R™ e e p tric - (The exceptional [geometric average] case might possibly 
simply be dismissed for serious consideration on the basis of numerical instabilities, with 
the associated volumes for this metric appearing to be actually infinite in nature. Also, 
as we will see below, Q KM is subject to particular severe jumps, perhaps decreasing the 
reliability of the estimates. The other five are rather close to 2 — but it is also somewhat 
intriguing that three of the estimated monotone metric ratios are quite close to one another 
(« 1.945), and therefore perhaps a common value unequal to 2.) This ratio-of-ratios can 
be easily be rewritten — as explicated in the Introduction — to take the form of a ratio of 
separability probabilities. That is, Q metnc is equivalently the ratio of the probability of sep- 
arability (p^( me *" c ' 6 l) f or a u qubit-qutrit states to the conditional probability of separability 
(Pg mefric ' 5 ') for those states on the (rank-five) boundary of the 3 5- dimensional convex set. 

An interesting conjecture now would be that this ratio (f2 met ™ c ) is equal to the integral 
value 2, independently of the (monotone or HS) metric used to measure the volumes and 
hyperareas. If, in fact, valid, then there is presumably a topological explanation |2jj| for 
this. (We were able to quite readily reject the speculation that this phenomenon might be 
in some way an artifact of our particular experimental design, in that we employ, as previ- 
ously discussed, only for simple computational convenience, a 34-dimensional subsequence 
of the 35-dimensional Tezuka-Faure sequence — rather than an ab initio independent 34- 
dimensional Tezuka-Faure sequence for the calculation of the hyperareas.) 

We must observe, however, that all the seven metrics specifically studied here induce the 
same (Haar) measure over 30 of the 35 variables — that is, the 30 Euler angles parameterizing 
the unitary matrices 0, ^8| — but not over the five independent eigenvalues of the 6x6 
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density matrix. Therefore, it is certainly valid to point out that we have not considered all 
types of possible metrics over the 35-dimensional space, but have restricted attention only 
to certain of those that are not inconsistent with quantum mechanical principles. (S. Stolz 
has pointed out, in a personal communication that, in general, one could modify a metric 
in the interior away from the boundary and outside the separable states, without affecting 
the metric on the separable states, thus changing R^+nonsep without changing R™f ric , but 



obviously then also altering the ratio of ratios (proportio proportionum 65j] ) Q metric . But 
presumably such a modification would lead, in our context, to the volume element of the 
so-modified metric not respecting Haar measure [cf. |Q, App. A]].) 

The topology of the (N 2 — l)-dimensional convex set of N x N density matrices has 
been laid out by Zyczkowski and Slomczynski |dtI . sec. 2.1]. The topological structure is 
expressible as 

[U(N)/T N ] x G N , (29) 

where the group of unitary matrices of size N is denoted by U (N) and the unit circle (one- 
dimensional torus ~ U(l)) by T, while Gn represents an (N — l)-dimensional asymmetric 
simplex. It would appear, however, that the set of separable states lacks such a product 
topological structure (thus, rendering integrations over the set — and hence the computation 
of corresponding volumes — quite problematical). 

In Fig. El is plotted the deviations from the conjectured integral value of 2 of the cumula- 
tive estimates of the ratio (Q GKS ) — given in Table ITXl — of the two hyperarea-to- volume 
ratios for the GKS monotone metric, the numerator [R S Qj^ umsep ) of Vt GKS being based on 
the entirety of qubit-qutrit states, and the denominator (Rq^s) being based on the bound- 
ary qubit-qutrit states only. (All the succeeding plots of deviations from the conjectured 
integral value of 2 will be drawn to the same scale.) In Figures E3 EH and El we show the 
corresponding plots based on the Bures, Hilbert-Schmidt and Kubo-Mori metrics, respec- 
tively. We note that the cumulative estimates in this last plot were relatively close to 2, 
before a sudden spike in the curve drove it upward. The values for the (quite unrelated) 
Bures and HS-metrics are rather close to 2, which is the main factor in our advancing the 
conjecture in question. 

It would, of course, be of interest to study comparable ratios involving 6x6 density 
matrices of generic rank less than 5. We did not originally incorporate these into our 
MATHEMATICA Tezuka-Faure calculations (in particular, since we did not anticipate the 
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FIG. 9: Deviations from the conjectured value of 2 of the cumulative estimates of Rgks, the 
ratio of hyp erarea-to- volume ratios for the Grosse-Krattenthaler-Slater (GKS or "quasi-Bures" ) 
monotone metric 
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FIG. 10: Deviations from the conjectured value of 2 of the cumulative estimates of Rsures, the 
ratio of hyp erarea-to- volume ratios for the Bures monotone metric 
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FIG. 11: Deviations from the conjectured value of 2 of the cumulative estimates of Rhs, the 
ratio of hyp erarea-to- volume ratios for the Hilbert-Schmidt metric. This plot is particularly flat in 
character 
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FIG. 12: Deviations from the conjectured value of 2 of the cumulative estimates of Rkm, the ratio 
of hyperarea-to-volume ratios for the KM (Kubo-Mori) monotone metric 



TABLE X: Counterparts for the qubit-qubit case N = 4 of the ratios of separability probabilities, 
based on 400 million points of a Tezuka-Faure sequence 



metric 



r 



metric 
sep+nonsep 



^metric 
sep 



metric 



Bures 
KM 
arith 
WY 
GKS 
geom 
HS 



12.1563 
0.506688 
2.19634 
2.93791 
6.03321 
4.02853 • 10~ 16 
51.9626 



6.58246 
0.348945 
1.2269 
1.73028 
3.3661 
7.1263 • 10~ 16 
25.9596 



1.84676 
1.45206 
1.79015 
1.69794 
1.79234 
0.565304 
2.00167 



apparent metric-independent phenomenon, we have observed here). In sec. IVI (TU below. we 
have, however, subsequently pursued such analyses. 



2. The N = 4 qubit-qubit case 

We adapted our MATHEMATICA routine used so far for the scenario A" = 6, so that it 
would yield analogous results for N — A. Based on 400 million points of a new independent 
15-dimensional Tezuka-Faure sequence, we obtained the results reported in Table |Xj (We 
now use the lower-case counterparts of the symbols R and Q to differentiate the N = 4 case 
from the A^ = 6 one.) Here, once more, the ratios-of-ratios ( Ci ; me * r * c ) tend to show rather 
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FIG. 13: Deviations from the possible true value of 2 of the cumulative estimates of the ratio of 
the rank-4 Hilbert-Schmidt separability probability to the rank-3 separability probability in the 
N = 4 qubit-qubit case. Note the greatly reduced scale of the y-axis. 

similar values, with the two exceptional cases again being the geometric average metric 
(which we suspect — like the maximal (Yuen-Lax) monotone metric, from which it is partially 
formed — simply gives infinite volumes and hyperareas) and the somewhat unstable KM 
monotone metric (which now gives an atypically low value!). We were somewhat surprised 
that the Hilbert-Schmidt metric again gives, as for iV = 6, a value quite close to 2. In 
Fig. El we show (on a comparatively very fine scale) the deviations from 2 of the cumulative 
estimates of the ratio of the Hilbert-Schmidt separability probability for the rank-4 states 
to that for the rank-3 states. 

However, it now seems fairly certain that if there is a true common value for UJ metr " lc 
across the metrics, then it is not an integral one (and thus possibly not a topological expla- 
nation) . The theoretical values predicted by (|2l)j) and (|2lj|) for r®^+ s nonsep and r B p +nonsep are 
16384/(429tt) « 12.1566 and 30^ « 51.9615, respectively. Also, consulting Table 6 of our 
earlier study [3], we find that using the conjectured and known values for the qubit-qubit 
case (N = 4) presented there gives us uj Bures = 8192/(14197r) w 1.83763 and a somewhat 
similar numerical value u arith = 408260608/73153125tt » 1.77646. 

Let us also indicate in passing that this new independent Tezuka-Faure sequence yields 
estimates that are quite close to previously known and conjectured values. For example, the 
ratios of the estimates of S^Q Bures ^ and g& Bures ) to their respective known values are 1.0001 
and .9999. Further, the ratios of the estimates of T,^{f wres ' and yj^ Snres ) ( our estimate being 
0.138119) to their respective conjectured values are 1.0001 and .99999. 
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FIG. 14: Deviations from the value 2 of the cumulative estimates of the ratio of the Hilbert-Schmidt 
probability of having a positive partial transpose for the 9x9 density matrices of rank 9 to the 
probability in the rank-8 case 

Let us take this opportunity to note that our analyses here indicate that the conjectures 
given in Table 6 of [3J] for the 14-dimensional hyperareas — denoted B s and B s+n there - 
appear to have been too large by a factor of 8. 

3. The N = 9 rank-9 and-8 cases 

K. Zyczkowski has indicated to us that he has an argument, if not yet fully rigorous, to 
the effect that the ratio of the probability of rank-iV states having positive partial transposes 
to the probability of such rank-(iV — 1) states should be 2 independently of N. Some early 
analyses of our — based on a so-far relatively short Tezuka-Faure sequence of 126 million 
points in the extraordinarily high (80)-dimensional hypercube — gave us a Hilbert-Schmidt 
rank-9/rank-8 probability ratio of 1.89125. (The analogous ratios for the monotone metrics 
were largely on the order of 0.15. In these same analyses we also — for our first time - 
implemented, as well, the computable cross-norm criterion for separability J68j], and found 
that many more density matrices could not be ruled out as possibly separable than with the 
[apparently much more discriminating] positive partial transposition criterion. The Hilbert- 
Schmidt probability ratio based on the cross-norm criterion was 0.223149.). In Fig. [T4"l we 
show the deviations from 2 of the cumulative estimates of the Hilbert-Schmidt rank-9/rank-8 
ratio based on the positivity of the partial transpose. However, this plot seems very unstable, 
so we must be quite cautious (pending a much more extended analysis) in its interpretation. 
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TABLE XI: Estimated ratios of both rank-6 and rank-5 qubit-qutrit separability probabilities to 
rank-4 separability probabilities 



metric 



rank-6/rank-4 ratio 



rank-5/rank-4 ratio 



Bures 
KM 
arith 
WY 
GKS 
geom 
HS 



20.9605 
12.2764 
17.4245 
15.5015 
18.3778 
1.30244 • 10~ 7 
33.9982 



10.7858 
5.4603 
8.52308 
7.96527 
9.44474 
0.00115724 
16.7652 



4- The N = 6 rank-4 and rank-3 cases 

The principal analyses above have been concerned with the full rank (rank-6) and rank-5 
6x6 density matrices. We adapted our MATHEMATICA procedure so that it would analyze 
the rank-4 and rank-3 cases, in a similar fashion. Now, we are dealing with 31-dimensional 
and 26-dimensional scenarios, vis a vis the original 35- and 34-dimensional ones. 

In a preliminary run, based on 35 million points of corresponding Tezuka-Faure sequences, 
not a sin gle rank-3 separable 6x6 density matrix was generated. (The general results of 
Lockhart 69j] — based on Sard's Theorem — tells us that the measures of rank-2 and rank-1 
6x6 separable density matrices must be zero, but not rank-3, as it appears we have observed 
[or near to it].) At that stage, we decided to concentrate further in our calculations on the 
rank-4 case alone. 

In Table IXT1 we report results based on one billion points of a (new/independent) 31- 
dimensional Tezuka-Faure sequence, coupled with our estimates obtained on the basis of 
our principal analysis, using the before-mentioned seven billion points. We note that for the 
Hilbert-Schmidt metric, 33.9982 (2 x 16.9991) is quite close to integral. In Fig. we show 
the cumulative estimates of the ratio from the value 34. (Of course, if the ratio of the rank-6 
HS separability probability to the rank-5 HS separability probability is exactly, in theory, 
equal to 2, and the rank-6/rank-4 separability probability exactly 34, then the rank-5/rank- 
4 ratio should be 17. Since it is based on greater numbers of sampled separable density 
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FIG. 15: Deviations from 34 of the cumulative estimates of the ratio of the rank-6 separability 
probability for the HS-metric to the rank-4 separability probability 

matrices, we suspect the sample estimate of the rank-6/rank-4 HS separability probability 
may perhaps be superior to the [less closely integral in value — that is, 16.7652] rank-5/rank- 
4 estimate.) 

Though the convergence to the predicted Hilbert-Schmidt volume was quite good 
(99.9654% of that given by the Zyczkowski-Sommers formula eq. (5.3)], for N = 6,n = 
2), we were initially rather disappointed/surprised that the predicted value of the Bures vol- 
ume was inaccurate by some 25%. This indicated to us that either the numerics were much 
more difficult for the Bures computation, or there was a possible error in our programming 
(which we were unable to locate) or even the possibility that something was incorrect with 
the specific Sommers-Zyczkowski formula 1^, eq. (4.19)] we were using, 

bN > n r((d n + i)M n )■ (30) 

This last possibility, in fact, proved to be the case, as we found that their formula (4.19) 
did not agree (for cases other than n — 0, 1, N — 1) with the more general formula (5.15) - 
reproduced above as (0) — and that using the correct formulation (5.15) (which we found 
also agrees with (4.18) of 17]) with (3 = 2, n = 2, N = 6, our numerical deviation was 
reduced from 25% to a more acceptable/less surprising .1%. A rectified version of their 
formula (4.19), 



q (2,Bures) _ -d, t ^ +1)/2 TjN _ n _ l j\(j + 2w)! 

' - v " ~ l r(K + l)/2)^=° [(j+n)!]" [6L) 
has since been posted by Sommers and Zyczkowski (quant-ph/0304041 v3) after this matter 
was brought to their attention. 
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D. Well-fitting formulas for the Bures and Hilbert-Schmidt separable volumes and 
hyperareas 



1. The Bures case 

Proceeding under the assumption of the validity of our conjecture above (regarding the 
integral value of 2 for Q metnc ), computational experimentation indicates that the Tezuka- 
Faure quasi-Monte Carlo separable Bures results can be quite well fitted by taking for the 
separable Bures 35-dimensional volume 

Z%/ ures) = « 1.03447 ■ 1(T 19 , (32) 

and for the separable (rank-five) 34-dimensional hyperarea 

^fures) = (2 43 . 3 . fc^)-! w L45449 . 1Q -« j (33) 

where by £ we denote volumes and hyperareas of separable states and 

CBures = V86429867T = Vvr • 2 ■ 11 • 19 ■ 23 ■ 29 ■ 31 « 5210.83. (34) 

The pooled sample estimates, T,flf ures ^ and £g 2 [ BMres \ as indicated in the caption to Ta- 
ble IVHl are 1.0142 • 10" 19 and 1.4707 • 10~ 18 . 

Then, we would have the Bures probability of separability of the (generically rank-six) 
qubit-qutrit states as 

p lBures,rank- 6 ] = ■ ■ U C Burea _ Q. 00142278, (35) 

and the Bures probability of separability of the generically rank-five qubit-qutrit states, 
exactly (by our integral conjecture) one-half of this. (However, there appears to be no 
obvious way in which the formulas immediately above extend the analogous ones in the qubit- 
qubit separable case |3j , which were hypothesized to involve the "silver mean" , OA g = 
Thus, it does not seem readily possible to use the results here to, in any way, support our 
earlier conjectures.) 

We have also devised another set of exact Bures formulas that fit our data roughly as 
well as (|52~|l and These are 

to Biirpt) 3 2 ■ 11 • 19 • 23 • 29 • 317T 1Q 
Zff ures) = ^-^ » 1.03497 • 10- 19 (36) 
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and for the separable (rank-five) 34-dimensional hyperarea 

sg Bures) = (2 43 • 5 7 )- 1 « 1.45519 • 1(T 18 . (37) 

2. The Hilbert- Schmidt case 

Additionally, we can achieve excellent fits to our Hilbert-Schmidt estimates by taking for 
the separable (rank-six) 3 5- dimensional volume 

Y,ff S) = (2 45 • 3 • 5 13 • 7v / 30) _1 « 2.02423 • 10~ 25 (38) 

(the sample estimate being 2.05328 • 10~ 25 ) and for the separable (rank-five) 34-dimensional 
hyperarea 

Eg H5) = (2 46 • 3 • 5 12 )- 1 « 1.94026 • lO" 23 ( 39) 

(the sample estimate being 1.94869 • 10~ 23 ). This gives us a Hilbert-Schmidt probability of 
separability of the generically rank-six states of 

\HSrank-® 3 10 ■ 7 4 ■ ll 3 ■ 13 2 ■ 17 2 ■ 19 ■ 23 ■ 29 ■ 31^5 . s 

p [Hb,rank 6j = V_ _ Q Q263115, (40) 

Zi ' 7T 

approximately 18.5 times the predicted Bures probability ()35p. (As we have noted, the Bures 
separability probability appears to be the greatest among the monotone metrics.) An upper 
bound of 0.166083 « (0.95) 35 on p ( j ffS >™»*-«*l i s g i ven i n Appendix G of Q. 

Our simple excellent Hilbert-Schmidt fits here led us to investigate whether the same 
could be achieved in the qubit-qubit (N = 4) case, using the same 400 million point Tezuka- 
Faure sequence employed in sec. IVI C 21 This, in fact, seemed definitely doable, by taking 



^(2,HS) 
J 4,0 



(3 3 5 7 v / 3) _1 w 2.73707 ■ 10~ 7 , (41) 



(the sample estimate being 2.73928 ■ 10 7 ) and 

Eff s) = (3 2 5 6 )" 1 w 7.11111 • 10" 6 (42) 

(the sample estimate being 7.11109 ■ 10~ 6 ). The exact Hilbert-Schmidt probability of sepa- 
rability of the generically rank-4 qubit-qubit states would then be 



Pf S ' rank ~ i] = ^ " 3 " 7 [: l l " 13V ^ « 0.242379. (43) 
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(A lower bound for Eg i of 2567r 7 /29805593211675 » 2.65834- 10~ 8 — that is, the volume 
of a 15- dimensional ball of radius ~ — appears obtainable from the results of [13[, although 
it is not fully clear to this reader whether the argument there applies to the two-qubit case 
[m — 2], since the exponent y — 1 appears.) Of course, one would now like to try to extend 
these Hilbert-Schmidt results to cases N > 6. Also, let us propose the formula 

^'"F^" 1 ' 02084 - 10 " (44) 

(the sample estimate [sec. I VI CM] being 1.04058 ■ 10~ 19 ). 



VII. DISCUSSION 

In the main numerical analysis of this study (sees. IV1 Al and IV1BJ) . we have directly 
estimated 28 quantities of interest — seven total volumes of the 35-dimensional space of 
qubit-qutrit states, seven 34-dimensional hyperareas of the boundary of those states, and the 
same quantities when restricted to the separable qubit-qutrit states. Of these 28 quantities, 

t 4- • a(2,B«res) c (2,Bures) a (2,HS) , a (2,HS)\ • l l c 

tour (that is, Dg q , >->6 i , og q an d ) were precisely known from previous 

analyses of Sommers and Zyczkowski [3, Hi- It is interesting to observe that the Tezuka- 
Faure quasi-Monte Carlo numerical integration procedure has, in all four of these cases, as 
shown in the corresponding table and figures (Table U and Figs, d El El EJ), slightly but 
consistently underestimated the known values — more so, it seems, with the 35-dimensional 
volumes, as opposed to the 34-dimensional hyperareas. (So, in statistical terminology, we 
appear to have biased estimators. The very same form of bias — in terms of the Bures metric 

- was observed in the precursor analysis [9| to this one, based on an independent, shorter 
Tezuka-Faure sequence. Randomizing deterministic algorithms — such as the Tezuka-Faure 

- can remove such bias |7o|.) This suggests that we might possibly improve the accuracy of 
the estimates of the 24 unknown quantities by scaling them in accordance with the magnitude 
of known underestimation. Also, we have in our several tables only reported the results at 
the (seven-billion point) end of the Tezuka-Faure procedure. We might also report results at 
intermediate stages at which the estimates of the 4 known quantities are closest to their true 
values, since estimates of the 24 unknown quantities might arguably also be most accurate 
at those stages. 

Of course, as we have done, taking the ratios of estimates of the volumes/hyperareas of 
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separable states to the estimates of the volumes/hyperareas of separable plus nonseparable 
states, one, in turn, obtains estimates of the probabilities of separability [l| for the various 
monotone metrics studied. (Scaling the estimated volumes and hyperareas by the corre- 
sponding estimates for the Bures metric, as we have done in certain of the tables above for 
numerical convenience and possible insightfulness, would be inappropriate in such a process.) 
Among the metrics studied, the Hilbert-Schmidt metric gives the largest qubit-qutrit proba- 
bility of separability (~ 0.0268283), while the Bures metric — the minimal monotone one - 
gives the (considerably smaller) largest separability probability (~ 0.00139925) among the 
monotone metrics studied (and presumably among all monotone metrics). The (Yuen-Lax) 
maximal monotone metric appears to give a null separability probability. 

In jj], we had attempted a somewhat similar quasi- Monte Carlo qubit-qutrit analysis (but 
restricted simply to the Bures metric) to that reported above, but based on many fewer points 
(70 million vs. the 7 billion so far used here) of a (Halton) sequence. At this stage, having 
made use of considerably increased computer power (and streamlined MATHEMATICA 
programming — in particular employing the Compile command, which enables the program 
to proceed under the condition that certain variables will enter a calculation only as machine 
numbers, and not as lists, algebraic objects or any other kind of expression), we must regard 
this earlier study as entirely superseded by the one here. (Our pooled estimate of the Bures 
volume of the separable qubit-qutrit systems here [Table VII] is 1.0142 ■ 10 -19 , while in j^, 
following our earlier work for N = 4 8], we formulated a conjecture (jj, eq. (5)]) — in 
which we can now have but very little confidence — that would give [converting from the 
SD metric to the Bures] a value of 2~ 35 - (2.19053- 10~ 9 ) « 6.37528- 10~ 20 .) We also anticipate 
revisiting — as in sec. IVI C 21 — the N = 4 (qubit-qubit) case J3J with our newly accelerated 
programming methods, in a similarly systematic manner. 

Perhaps, in the future, subject to research priorities, we will add to the 7 billion points of 
the Tezuka-Faure sequence employed above, and hope to report considerably more accurate 
results in the future (based on which, possibly, we can further appraise the hypotheses 
offered above as to the values of the various volumes and hyperareas). Also, we may seek 
to estimate the hyperarea of that part of the boundary of the separable qubit-qutrit states 



consisting of generically rank-six 6x6 density matrices [3j, |63J, though this involves a much 
greater amount of computation per point. (This would entail first finding the values, if any, 
of the undetermined [35-i/i] parameter that would set the determinants of the two forms of 
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partial transpose equal to zero, and then — using these values — ascertaining whether or 
not all the six eigenvalues of the resultant partial transposes were nonnegative.) 

In this study, we have utilized additional computer power recently available to us, to- 

S'her with an advanced quasi- Monte Carlo procedure (scrambled Faure-Tezuka sequences 
- the use of which was recommended to us by G. Okten, who provided a cor- 
responding MATHEMATICA code). Faure and Tezuka were guided "by the construction 
C (i) = A d)p(i-i) and by 

some possible extensions of the generator formal series in the frame- 
work of Neiderreiter". (A^ % > is an arbitrary nonsingular lower triangular [NLT] matrix, P 



is the Pascal matrix 



Ilj and C w is a generator matrix of a sequence X.) Their idea was 



to multiply from the right by nonsingular upper triangular (NUT) random matrices and 
get the new generator matrices = p( l-1 )[/W for (0, s)-sequences (45! I^. "Faure-Tezuka 
scrambling scrambles the digits of i before multiplying by the generator matrices . . . The 
effect of the Faure- Tezuka-scrambling can be thought of as reordering the original sequence, 
rather than permuting its digits like the Owen scrambling . . . Scrambled sequences often have 
smaller discrepancies than their nonscrambled counterparts. Moreover, random scramblings 
facilitate error estimation" 2ij, p. 107]. 

It would be interesting to conduct analogous investigations to those reported here (N = 6) 
and in [^] for the case N = 4, using quasi-random sequences other than Tezuka- Faure ones 
particularly those for which it is possible to do statistical testing on the results (such 
as constructing confidence intervals) 170]. It is, of course, possible to conduct statistical 
testing using simple Monte Carlo methods, but their convergence is much weaker than that 
of the quasi- Monte Carlo procedures. Since we have been dealing with extraordinarily high- 
dimensional spaces, good convergence has been a dominant consideration in the selection of 
numerical integration methodologies to employ. 

"It is easier to estimate the error of Monte Carlo methods because one can perform a 
number of replications and compute the variance. Clever randomizations of quasi-Monte 
Carlo methods combine higher accuracy with practical error estimates" 0, p. 95]. G. 
Okten is presently developing a new MATHEMATICA version of scrambled Faure-Tezuka 
sequences in which there will be a random generating matrix for each dimension — rather 
than one for all the dimensions together — which will then be susceptible to statistical 
testing |to| . 

At the strong urging of K. Zyczkowski, we disaggregated the pooled results in 
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the last column of Table IX into the part based on partial transposition of four 
three-by-three blocks and obtained {1.966,2.53505,2.04826,1.94679,1.96481,9.32089 x 
10~7, 1.99954} and into the part based on nine two-by-two blocks and obtained 
{1.91846, 1.91976, 2.04, 1.94539, 1.92476, 0.0001227, 2.05803}. We bring the attention of the 
reader to the particular closeness to 2 of the first (Hilbert-Schmidt) ratio. 
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